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I. INTRODUCTION 



The standard way to analyze the high-energy behavior of pQCD amplitudes is the direct summation of Feynman 
diagrams. In the Regge limit s » t^m? (where to is a characteristic mass or virtuality of colliding particles) there 
is an extra parameter in addition to as , namely as In ^ . At pre-asymptotic energies we may have a window where 
the leading logarithmic approximation (LLA) 

g 

as < 1, as In ^ 1 (1) 

TO 
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is valid. Following the pioneering work by Lipatov the amplitude in this region of energies was shown to be 
determined by the BFKL pomeron [2] which leads to the behavior of QCD cross sections of the type 

-i'^) = du f-i F{.) ^ f— (2) 



where is a number of colors (iVj, — 3 for QCD), uj{v) — ^Nc[2il){\) — — tp{^— iv)] is the BFKL intercept, 

and F{i>) is the "pomeron residue" which depends on the process. 

The power behavior of the cross section ([2]) violates the Froissart bound a < In^ s, and therefore the BFKL pomeron 
describes only the pre-asymptotic behavior at intermediate energies when the cross sections are small in comparison to 
the geometric cross section 2ttR^. In order to find the exact window of the BFKL applicability one needs to calculate 
the first correction to the BFKL amplitude. 

The importance of corrections to LLA result in high-energy QCD is twofold. As I mentioned above, to get the region 
of application of a leading-order result one needs to find the next-to-leading order (NLO) corrections. However, in the 
case of the BFKL amplitude ^ there is another reason why NLO corrections are essential. Unlike for example the 
DGLAP evolution, the argument of the coupling constant in Eq. © is left undetermined in the LLA, and usually it is 
set by hand to be of order of characteristic transverse momenta. Careful analysis of this argument is very important 
from both theoretical and experimental points of view, and the starting point of the analysis of the argument of as 
in Eq. ^ is the calculation of the NLO BFKL correction. 

The NLO correction to the BFKL pomeron intercept uj{i^) was found by two groups - Lipatov and Fadin f§\ and 
Ciafaloni and Camichi [4] - after almost ten years of calculations. The result has the form 



TT 



An 



b 



2Nr 



67 TT^ lOUf 



Si'') = - TTTFx'i^) + 6C(3) + TT - ^ - 77^ + x"(^) - 2$ (^.) - 2<i>(-^.) (3) 
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where b = -yiVc — xi'^) = ^ 2C — iJj{^ + iv) — — iv) (C = — ^/'(l) is Euler's constant), and 

$H = -^'-^t-H-[!^ + 2Li2W]. (4) 

The term proportional to b depends on the choice of the argument of coupling constant. The choice of this term in 
the r.h.s. of Eq. ([3]) corresponds to a{qiq2) where qi are momenta of scattered BFKL gluons, see the discussion in 
Ref. . The pomeron intercept is related to anomalous dimensions of the twist-2 gluon operators - it determines the 
asymptotics of anomalous dimensions 7j as — 1 = w — )■ 0. It is worth noting that the result of explicit calculation 
of 3-loop anomalous dimensions at the 3- loop level Q agrees with Eq. ([3]). As we see from Eq. ([3]) the full NLO 
description of the amplitude implies the knowledge of the "pomeron residue" F{v) at the NLO level. A classical 
example is the scattering of virtual photons in QCD where F{v) is proportional to product of two amplitudes of 
emission of two t-channel gluons by the upper and the lower virtual photon. In the leading order, these "impact 
factors" are known for a long time since they differ from the corresponding QED result only by trivial color factors. 
However, to the best of my knowledge, there is no complete analytical calculation of the NLO photon impact factor 
in the literature (for a combination of numerical and analytical results, see Ref. @). I will present the analytical 
result for the photon impact factor in the coordinate space which is compact and conformally invariant; however, 
the Fourier transformation to the momentum space is not performed yet. The calculations use different approach to 
high-energy scattering based on the operator expansion in Wilson lines Q rather than on the direct summation of 
Feynman diagrams. 

A general feature of high-energy scattering is that a fast particle moves along its straight-line classical trajectory 
and the only quantum effect is the eikonal phase factor acquired along this propagation path. In QCD, for the fast 
quark or gluon scattering off some target, this eikonal phase factor is a Wilson line - the infinite gauge link ordered 
along the straight line collinear to particle's velocity n'^: 



{x±) ~ Pexp-^ — ig I du Tip A^[un + x±) >, (5) 

J — oo 



Here is the gluon field of the target, xa_ is the transverse position of the particle which remains unchanged 
throughout the collision, and the index Y labels the rapidity of the particle. Repeating the above argument for the 
target (moving fast in the spectator's frame) we see that particles with very different rapidities perceive each other as 
Wilson lines and therefore these Wilson-line operators form the convenient effective degrees of freedom in high-energy 
QCD (for a review, see ref. [1]). 
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As an example of high-energy process let us consider the deep inelastic scattering from a hadron at small Xb = 
Q"^ I {2p ■ q). The virtual photon decomposes into a pair of fast quarks moving along straight lines separated by some 
transverse distance. The propagation of this quark-antiquark pair reduces to the "propagator of the color dipole" 
U{x±)W {y±) - two Wilson lines ordered along the direction collinear to quarks' velocity. The structure function of a 
hadron is proportional to a matrix element of this color dipole operator 

ll''{x^,y^) = 1 - l-Tr{U^{x^)U^^{y^)} (6) 

switched between the target's states {Nc — 3 for QCD). The gluon parton density is approximately 

XBGixB,^l^ ^Q^) ^ {p\U''{x^,0)\p) (7) 

where Y — In^. (As usual, we denote operators by "hat"). The energy dependence of the structure function is 
translated then into the dependence of the color dipole on on the rapidity Y . There are two ways to restrict the 
rapidity of Wilson lines: one can consider Wilson lines with the support line collinear to the velocity of the fast- 
moving particle or one can take the light-like Wilson line and cut the rapidity integrals "by hand" (see Eq. ([50)1 
below). While the former method appears to be more natural, it is technically simpler to get the final results with 
the latter method of "rigid cutoff" in the longitudinal direction. 

Eq. ([71 means that the small-x behavior of the structure functions is governed by the rapidity evolution of color 
dipoles [3, [13] • At relatively high energies and for sufficiently small dipoles we can use the leading logarithmic 
approximation ([T]) [ug aslnxs^l) and get the non- linear BK evolution equation for the color dipoles [tI. [Tl|: 

-^Z^^(zi,Z2) = ^ fd^z,^^[U''{z,,z,)+U''{z,,z,))-U''{z,,zs)~U^iz,,Z3)U''{z,,z,)] (8) 
ai zvr J ^13^23 

where 1" = In ^ and zi2 = zi — Z2 etc. The first three terms correspond to the linear BFKL evolution and describe 
the parton emission while the last term is responsible for the parton annihilation. For sufficiently low xb the parton 
emission balances the parton annihilation so the partons reach the state of saturation [12] with the characteristic 
transverse momentum Qs growing with energy 1/xb (for reviews, see Ref. [13]) 

It is easy to see that the BK equation ([8]) is conformally invariant in the two-dimensional space. This follows from 
the conformal invariance of the light-like Wilson lines. The Wilson line 

U{x_i_) = Pexp {-ig f dx+ A+{x+,x±)\ (9) 

J — oo 

is invariant under the inversion x^ — >■ x^ jx^ (with respect to the point with zero (-) component). Indeed, (a;+, .tj^)^ = 
x\ so after the inversion x^ — > x\^lx\ and — > x'^ lx\ and therefore 



V{xx_) ^ Pexp -^5/ d— A+{ — ,xx_)\ = U{x^/x\) (10) 

Moreover, it is easy to check formally that the Wilson line operators lie in the standard representation of the conformal 
Mobius group SL(2,C) with conformal spin (see |l3l))- It should be noted that the conformal invariance of the linear 
BFKL equation was first proved in Ref. [T5| . 

The NLO evolution of color dipole in QCD is not expected to be Mobius invariant due to the conformal anomaly 
leading to dimensional transmutation and running coupling constant. To understand the relation between the high- 
energy behavior of amplitudes and Mobius invariance of Wilson lines, it is instructive to consider the conformally 
invariant M = A super Yang-Mils theory. This theory was intensively studied in recent years due to the fact that 
at large coupling constants it is dual to the IIB string theory in the AdSs background. In the light-cone limit, the 
contribution of scalar operators to Maldacena- Wilson line [16] vanishes so one has the usual Wilson line constructed 
from gauge fields and therefore the LLA evolution equation for color dipoles in the N — A SYM has the same form 
as At the NLO level, the contributions from gluino and scalar loops enter the game. 

As I mentioned above, formally the light-like Wilson lines are Mobius invariant. However, the light-like Wilson 
lines are divergent in the longitudinal direction and moreover, it is exactly the evolution equation with respect to this 
longitudinal cutoff which governs the high-energy behavior of amplitudes. At present, it is not known how to find 
the conformally invariant cutoff in the longitudinal direction. When we use the non-invariant cutoff we expect, as 
usual, the invariance to hold in the leading order but to be violated in higher orders in perturbation theory. In our 
calculation we restrict the longitudinal momentum of the gluons composing Wilson lines, and with this non-invariant 
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cutoff the NLO evolution equation in QCD has extra non-conformal parts not related to the running of coupling 
constant. Similarly, there will be non-conformal parts coming from the longitudinal cutoff of Wilson lines in the 
W = 4 SYM equation. I will demonstrate below that it is possible to construct the "composite conformal dipole 
operator" (order by order in perturbation theory) which mimics the conformal cutoff in the longitudinal direction so 
the corresponding evolution equation is Mobius invariant. With the NLO accuracy this composite operator has the 
form 14] (a is an arbitrary constant): 

mu^M}]l°;' (11) 

ZTT J ^13-^23 ''^13-^23 

where the second term is a "counterterm" restoring the conformal invariance lost because of the cutoff (pO)) . This is 
quite similar to the construction of the composite renormalized local operator in the case when the UV cutoff does not 
respect the symmetries of the bare operator - in this case the symmetry of the UV-regularized operator is preserved 
order by order in perturbation theory by subtraction of the symmetry-restoring counterterms. 

The NLO BK equation written in terms of these composite conformal operators is conformally invariant (in A/" = 4 
SYM) and so are the impact factors [14]. Below I present the results for these impact factors and for the "pomcron 
residue" F{u) in a simple case of correlator of four scalar currents in the Regge limit. 

In QCD wc do not have the conformal invariance so the corresponding composite operators are not, strictly speaking, 
conformal. Still, if we write down the NLO BK equation in terms of these operators it has a nice property of being a 
sum of conformal part and the running-coupling part proportional to 6 = — ^Uf. I calculate the NLO coefficient 
function in the expansion of two electromagnetic currents in these composite operators which determines the photon 
impact factor at the NLO level. 

The paper is organized as follows. In Sect. II I remind the general logic of the operator expansion using the example 
of light-ray OPE. In Sect. Ill I formulate the program of high-energy OPE and carry it out in subsequent Sections. 
In Sect IV we calculate the high-energy amplitudes in A/" = 4 SYM in the leading order and then in the NLO. Sect 
V is devoted to the NLO high-energy amplitudes in QCD and Sect. VI contains the conclusions. The details of 
the calculation of Feynman diagrams for the NLO evolution of color dipoles are not presented here (see the original 
publications [ij, [l7l - [l9| ) but we will outline the calculation of impact factors in Af — A SYM and in QCD. 



II. THE LOGIC OF OPE 



Let me first remind the logic of usual operator expansion near the light cone. A typical example is the calculation of 
the structure functions of deep inelastic scattering (DIS) at moderate xb — Q^/{2p- q) determined by the T-product 
of two electromagnetic currents T{j^{x)jy{y)} switched between the target states. The first step is to identify the 
relevant operators. To this end, we formally put to infinity, then {x — y)^ = and we see that the relevant 
operators are the light-ray ones of the type of ip{x)^^[x, y]'il'{y) where 

[x,y] = Pex^p^-igJ du {x - y)'' Af,{ux - uy + y)^ (12) 

is a standard notation for a straight-line gauge link connecting points x and y. 

Now, to get the behavior of structure functions one needs to perform the following four steps: 

• Separate relevant Feynman loop integrals over k± in two parts - the coefficient functions (with transverse 
momenta fc^ greater than the factorization scale /i^) and parton densities - matrix elements of light-ray operators 
with fc^ < /i^ (see Fig. 1). Technically, the cutoff of the transverse momenta in these matrix elements is done 
by adding counterterms with normalization point /i^ . 

• Find the evolution equations of light-ray operators with respect to /x^. 

• Solve of the corresponding evolution equation(s). 

• Assemble the result for structure functions: take some initial conditions at low normalization point (usually 
around /x^=lGeV^), evolve parton densities to high and multiply the result by the coefficient function. 

Let me explain these steps at the NLO level in detail. 
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FIG. 1: Operator product expansion near the light cone. Gauge hnk is denoted by a dotted hne. 



One can calculate the coefficient functions by direct computation of Feynman diagrams with quark and gluon tails 
or by using the Feynman diagrams in external quark and gluon fields (see e.g. [l^, [HI)- A typical NLO contribution 
coming from the upper part of last diagram in Fig. [T] has the form 





1- 


/ du 


Jo 





where 



is a standard "plus" prescription 



[ipix)[x, ux] ytipiux)]^-^- (l + ^[InxV^ + In 



(13) 





1- 




/ du 


fin) ^ f 


/o 


-U. 


+ Jo 



(14) 



and the leading-twist light-ray operator [^{x)[x,ux] ^i/j{ux)Y'*^' is taken at x'^ = 0. Formally, one can consider 
the operator at x where x is close to x and light-like. Alternatively, one can define a "leading twist" prescription by 
subtraction of all higher twists in the non-local form, see In any case, the light-ray operator [^{x)[x, ux] ;j6^{ux)Y'*" 
has specific light-cone UV divergencies (in addition to usual UV contributions assembled to the running coupling 
constant) which are regularized by adding light-ray counterterms. The resulting composite operator (original operator 
minus counterterms) depends on the renormalization point /i and the renorm-group equation for these operators leads 
to the DGLAP evolution equation for parton densities. 

A typical contribution to the NLO evolution equation of the quark light-ray operator coming from the bottom part 
of diagram in Fig. [T] has the form 



f,—[^p{x)[x,o] MOW,. 



ajji) 
2n 



-Cp 





1- 


1 du 


Jo 





[tp{x)[x,ux] ;/:^{uxW^ (l + — [Inx^^^ + In u] 
■ V 47r 



(15) 



{cp = {N^ — l)/2Nc). This completes the second step of the above four-step program. 

The third step is the solution of the evolution equation of (jlSp type. I will present this solution in a simple case of 
forward matrix elements of quark operators for a non-singlet case determined by quark parton densities. For example, 
in the case of unpolarizcd proton we get 



{pmx)[x,o]Mo)-qAO:^]qi^)tt.\p) = 2{px) 



) / duj Vq{u)[e'P-'"^ 
Jq 



-tp-xujl 



(16) 



where p is the proton momentum and Dq{uj) is a parton density of quark q (u, d, or s). The solution of the evolution 
equation for parton densities is given by the Mellin integral. In terms of light-ray operators it has the form 



(^(x)[x,o] i^x'^mtt.) 



— ^ — 'too 



/_a^(^\ — f- 
\a^{uo) / 



du u'{il^{ux)[ux,G\ ^A"V(0)]rt.> 

(17) 

where A" is a flavor Gell-Mann matrix and (O) means any forward matrix element of the operator O. Also, 7j(as) — 

if^j^ + 7^^'' + ••■ is the anomalous dimension of the light-ray operator J^du u^-)jj{ux)[ux, 0] /A''-0(O) and bi is 
the second coefficient of Gell-Mann-Low function. 
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Now, to perforin step four and get a structure function {—3Fi + in our case) one can take parton densities at 
low /iQ ~ IGev^, put it into Eq. (IT7|) . take fi^ = and multiply the result by the coefficient function (|13p taken at 
fi^ = Q^. In future, one should get the parton densities at low fj? from non-perturbative (lattice) QCD but at present 
people use well established models. The evolution of parton densities in QCD is calculated up to the third order in 
as [5| and the results agree with experiment in a broad range of parameters and xb ■ 

III. HIGH-ENERGY OPE IN WILSON LINES 

Now we want to extend the four-step program from the previous section to describe the high-energy amplitudes. A 
typical process is deep inelastic scattering at small values of Bjorken x. Since we are interested now in xb evolution 
rather than evolution, as step one we factorize in rapidity instead of transverse momenta. As a preliminary step we 
identify the relevant operators. The virtual photon splits into quark-antiquark pair, and if we set formally the energy 
of incoming photon to infinity we see that these quark and antiquark travel along the light-like classical trajectories. 
This is the well-known general result from quantum mechanics - the fast particle moves along its straight-line classical 
trajectory and the only quantum effect is the eikonal phase factor acquired along this propagation path. In QCD, for 
fast quark or gluon scattering off some target, this eikonal phase factor is a Wilson line ([5]) - an infinite gauge link 
ordered along the straight line collinear to particle's velocity n'^: 




FIG. 2: High-energy operator expansion in Wilson lines 



Having identified the operators, we can repeat the 4 steps which give us the high-energy (small- a; for DIS) evolution 
of the amplitude. 

• Separate the relevant Feynman loop integrals over the longitudinal momentum a in two parts - the coefficient 
functions ( "impact factors" ) with a greater than the rapidity factorization scale a = and matrix elements of 
Wilson-line operators with a < a (see Fig. We were not able to find the analog of dimensional regularization 
for longitudinal divergence so we cut the integration over a < a "by hand" . 

• Find the evolution equations of color dipoles with respect to the cutoff in rapidity Y . 

• Solve of the corresponding evolution equation(s). 

• Assemble the result for structure functions: take the initial conditions at low energy, evolve color dipoles to high 
rapidity Y ~ Ya and multiply the result by the corresponding impact factor. 

Unlike the DGLAP evolution discussed above, the BK evolution equation for color dipoles with respect to energy 
is non-linear which makes it considerably more difficult to solve. At present, for the experimentally interesting case 
of DIS from proton or nucleus we have no analytical solution so one has to rely upon the approximate solutions and 
numerical simulations. One can linearize the BK equation and get the usual linear BFKL evolution but the validity 
of this linearization in the saturation region of small x is questionable. Still, there are purely perturbative processes 
where the BFKL equation gives the correct (pre-asymptotic) behavior at high energies: for example the scattering of 
virtual photons with equal (and high) virtualities, or observation of two Mueller-Navelet jets [23| in hadron-hadron 
scattering. In this case, the leading-order analysis has been performed, but the generalizations of the LO results to 
the next-to-leading order has not been obtained since it is difficult to take into account the running-coupling effects 
in the BFKL equation. 
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By the same token, since in = 4 SYM the couphng does not run, the NLO BFKL program can be performed 
to the very end. I will present the NLO analysis of high-energy "scattering of two scalar currents" and get the NLO 
result for this correlation function in an explicit form. In subsequent two chapters we will carry out steps (1-4) of our 
program for A/" = 4 SYM and then for QCD (where we'll discuss only steps 1 and 2). 



IV. HIGH-ENERGY AMPLITUDES IN A/" = 4 SYM IN THE NEXT-TO-LEADING ORDER 

A. Regge limit and Pomeron in — 4 SYM 

As we mentioned above, at first we will find the NLO amplitudes at high energies for A/" = 4 SYM and turn to 
QCD later. We use the A/" = 4 Lagrangian in the form (see e.g. Ref. [23|): 



1 



2 fabc fhnc^a^b^^l^^m 

^sAB^s\a rate 



(18) 



Here $j are scalars, A"'^ gluinos and EJj = {YXb ' ''■'^ab) ^ ^ij = O^abt ^^^ab) where are standard 't Hooft 
symbols. The bare propagators are 



{\f{x)\'^{y)) 



(19) 



and the vertex of gluon emission in the momentum space is proportional to {ki — k2)^T'^5jj for the scalars and a'^T"' 
for gluinos. Here cr^ — (1, a), = (—1, (?) where a are usual Pauli matrices and our metric is g^" — (—1, 1, 1, 1). 
For simplicity, let us consider correlation function of four scalar currents 



A{x,y,x',y') = {x - y)\x' - y'f{0{x)OHy)0{x')OHy')) 
where O = ^^^^Tr{Z^} {Z = "^('^'i + ^^2)) is a renorm-invariant chiral primary operator. 



(20) 



In a conformal theory this four-point amplitude A{x, y; x', y') depends on two conformal ratios which can be chosen 



as 



f\2 



R 



{x-x'f{y-y') ^ 
[x - yY{x' - y'Y' 



R 



{x~y'f{v-x'? 1 
{x~x'Y{y-y'Y R 



(21) 



We are interested in the behavior of the correlator pO)) in the high-energy (Regge) limit. In the coordinate space it 
can be achieved as follows: 



2 2 

x = px^-pi+x^, y^py*-pi+y^, 



P X, -P2 

s 



y 



p'y',-P2 + y± 



(22) 



with p, p' 00 and a;* > > y*, a;^ > > y^. (Strictly speaking, p 00 oi p' 00 would be sufficient to reach the 
Regge limit). Hereafter I use the notations x, = —p^x^, = —p2X^i where pi and p2 are light-like vectors normalized 
by —2{pi,p2) — s. These "Sudakov variables" are related to the usual light-cone coordinates x^ — -^(a;" ± x^) by 

X* — x^y^s/2, X, ~ x^ so X — |x*pi -I- |x,p2 + x±. The metric is sr''''=(-l, 1,1,1) so x^ = —|x,x* -|- x^^. In the 

Regge limit ([22l) the full conformal group reduces to Mobius subgroup SL(2,C) leaving the transverse plane (0, 0, z±) 
invariant. 

As demonstrated in Ref. ^25j . the pomeron contribution in a conformal theory can be represented as an integral 
over one real variable v 



(x-yrix'-yr{Oix)OHy)0{x')OHy')) - 7^ di^ h{^) 



tanh'; 



-F(i/)0(r,i/)i?^^('') 



(23) 



Here oj{v) = uj{0,v) is the pomeron intercept, f+{v) = /+(w(z^)) where f+{uj) = (e™'^ ~ 1)/ siuTrw is the signature 
factor in the coordinate space, and F(v) is the "pomeron residue" (strictly speaking, the product of two pomeron 
residues). The conformal function i7(r, i^) is given by a hypergeometric function (see Ref. [26| ) but for our purposes 
it is convenient to use the representation in terms of the two-dimensional integral 



17(r, v) ^ ^\ <rz 



{-2k -Cf 



{-2k'-Q^ 



(24) 
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where ( = ^ + z'\p2 + z±^ and 



, \/s ,Pl ,2 , , VS /Pi , ,2 , . 

are two SL(2,C)-invariant vectors [1^ (see also [l^]) such that — '^'^ — 



(25) 



4a:, a. 



•^^T , ^, and therefore 
4a;'. yl 



2 /2 



1 



4(k-k')' 



r 



(26) 



In our limit x — x^, x' = x' j_ and similarly for y. Note that all the dependence on large energy (= large p, p') 
is contained in R^^^'^\ 

The dynamical information about the conformal theory is encoded in two functions: pomeron intercept and pomeron 
residue. The pomeron intercept is known both in the small and large as limit. The NLO intercept at small Us was 
calculated in A/" = 4 SYM by Lipatov and Kotikov [111 



M = ^^c(x(^) + ^[6C(3)-yx(^)+x"(^)- 2$(^.)-2$(-;.)]) 



(27) 



Our main goal is the description of the amplitude in the next-to-leading order in perturbation theory, but it is worth 
noting that the pomeron intercept is known also in the limit of large 't Hooft coupling A — AiragNc 



2-2- 



(28) 



where 2 is the graviton spin and the first correction was calculated in Ref . |29l. [3 01 . 

The pomeron residue F{v) is known in the leading order both at small [2^3^321 and large [1^ 't Hooft coupling 



n/ \ >-+0 ,2 TI'SiuTTl' A^oo TT^ ( 1 -f J^) ^ 

^(^) ^ ^ A 3 ' -^(^) ^ 



sinh -Kv 



(29) 



To find the NLO amplitude, we must also calculate the "pomeron residue" F{v) in the next-to-leading order. In the 
rest of this section we will do this using the four steps of the high-energy operator product expansion in Wilson lines. 



B. High-energy OPE in the leading order 



1 . Leading order: impact factor 



As I discussed above, the main idea behind the high-energy operator expansion is the rapidity factorization. At 
the first step, we integrate over gluons with rapidities Y > i] and leave the integration over Y < for later time, see 
Fig. 2. It is convenient to use the background field formalism: we integrate over gluons with a > a and leave gluons 
with a < cr as a background field, to be integrated over later. The result of the integration is the coefficient function 
("impact factor") in front of the Wilson-line operators with rapidities up to 77 = In cr: 



U;^ = Pexp 



^9 



\ak\)e'^--A^{k) 



(30) 



where the Sudakov variable ak is defined as usual, k = akPi + PkP2 + ■ The impact factor is given then by a set of 
Feynman diagrams in the external field of gluons with a < a. Since the rapidities of the background gluons are very 
different from the rapidities of gluons in our Feynman diagrams, the background field can be taken in a form of the 
shock wave due to the Lorentz contraction. It is very easy to derive the expression of a quark (or gluon) propagator 
in this shock-wave background. We represent the propagator as a path integral over various trajectories, each of them 
weighed with the gauge factor Pexp(i(7 Jdx^A^^) ordered along the propagation path. Now, since the shock wave is 
very thin, quarks (or gluons) do not have time to deviate in transverse direction so their trajectory inside the shock 
wave can be approximated by a segment of the straight line. Morover, since there is no external field outside the 
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FIG. 3: Propagator in a shock- wave background 

shock wave [i^ the integral over the segment of straight line can be formally extended to ico limits yielding the 
Wilson-line operator. 

Thus, the structure of the propagator in a shock-wave background is as follows: 
[Free propagation from initial point x to the point of intersection with the shock wave z] 
X [interaction with the shock wave described by the Wilson-line operator Uz] 
X [Free propagation from point of interaction z to the final point ?/] . 



The explicit form can be taken from Ref. JlJ 



1 



di"'' 

47r2 [{x - zY + ie] '-^ At:^ [{z - yY + ie] 



2XU 



s'6 
64:Tr^x*y* Jq 



(31) 



where Z = — -^(k • Q = —^{x — y), + ^^''-^ — '•^ ^^-'^ . Note that the interaction with the shock wave does not 
change the a-component of the scalar particle's momentum. 

Let us calculate the impact factor taking x, = y, = for simplicity. The leading-order impact factor is proportional 
to the product of two propagators ((3T|) . see Fig. IH 




FIG. 4: Impact factor in the leading order. 



where 



7^ = 



(a; - yf 



x^y^ZiZ2 



(32) 



(33) 



and C,i = ^ + zfp2 + 2,^ , Zi = —^{k-Q) — . Note that the leading-order impact factor is conformally 

(Mobius) invariant - it goes into itself under the inversion ([TO)) . 
This formula can be promoted to the operator equation as follows 



(^x-yfT{0{x)0{y)}^^ = 



1 



7r2(7V2 _ 1) 



'12 



(34) 



where u a ~ ^^l^^j'J is the characteristic a's in the scalar loop which serve as an upper bound for rapidity of Wilson- line 
gluons. (Recall that and are of the same order of magnitude as seen from Eq. ^2f2\ ). 
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We'll need later the proiection of the T-product in the l.h.s. of this equation onto the conformal eigenfunctions of 
the BFKL equation [l5| 



Eu,n{ziO, Z20) 



' Z12 ■ 




■ Z12 ' 






■Z10Z20- 



(35) 



(here z = Zx-\- izy, z = Zx — izy, zio = Zi — zq etc.). Since Cs are scalar operators, the only non-vanishing contribution 
comes from projection on the eigenfunctions with spin 0: 



1 f dzidz2_ ^2 



Z12 



'12 



^2 2 



(-2k • Co)^ J r(l - 2w) cosh TTv 



(36) 



where Co = ^ + zl^p2 + zq^. 

Now, using the decomposition of the product of the transverse (5-functions in conformal 3-point functions 

E^^„{ziq,Z2q) [11 



5^^\z,-u^,)5^-\z2-W2) ^ E r^'^^ d^PK, 

„ — c^J ^ ^12"'l2 J 



{Wi - p,W2 - p)E^^riizi - P,Z2 - p) (37) 



we obtain (dots stand for contributions of higher spins n which we do not need for our correlator (j2Qp ) 



ix-y)'T{d{x)&iy)} 



dv 1 d zq 



47rcosh7ri/ r(l - 2w) \{-2k-C,oY 



(^^2^)^">^(.,.o) + ... (38) 



where 



- 1 fd zid Z2 

U (i/,2:n) = — 



{-^y^"' w{z,,z2) 



(39) 



and W^{zi, Z2) is a "color dipole in the adjoint representation" 



U^z^,Z2) = I - j^T^Tr{U:^Ul:} 



(40) 



2. Leading order: BK equation 

Next step is to obtain the evolution equation for color dipoles in the leading order in as- For the light-cone 
dipoles, the contribution of scalar operators to Maldacena- Wilson line [1^1 vanishes so one has the usual Wilson line 
constructed from gauge fields and therefore the LLA evolution equation for color dipoles in the N = A SYM has the 
same form as in QCD. 

To find the evolution of the color dipole © with respect to rapidity of the Wilson lines in the leading log approx- 
imation we consider the matrix element of the color dipole between (arbitrary) target states and integrate over the 
gluons with rapidities Yi > F > I2 = i^i — leaving the gluons with F < F2 as a background field (to be integrated 
over later) . In the frame of gluons with Y ^ Yi the fields with Y <Y2 shrink to a pancake and we obtain the four 
diagrams shown in Fig. [S] Technically, to find the kernel in the leading-ordrer approximation we write down the 
general form of the operator equation for the evolution of the color dipole 

^Triillul-^} ^ K^o^r{Ulull} + ... (41) 

(where dots stand for the higher orders of the expansion) and calculate the l.h.s. of Eq. (|4T|) in the shock-wave 
background 

^(Tr{J7];;7lf }),hockwavc = (ifLoTr{?7];;7l^^})shockwavo (42) 
In what follows we replace (...) Shockwave by (...) for brevity. 
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FIG. 5: Leading-order diagrams for the small- a: evolution of color dipole. Gauge links are denoted by dotted lines. 

With future NLO computation in view, we will perform the leading-order calculation in the lightcone gauge P2 ~ 
0. The gluon propagator in a shock-wave external field has the formflTl. Issj 

2 J TV^yix — zy + ie\'^ dl ■k'^[{z — yy + ie\'^ 

where ^ can be either g or g-z^ which leads to the same result. (This is obvious for the leading order and 
correct in NLO after subtraction of the leading-order contribution, see Eq. (|67|) below). 
We obtain 

poo fO i*e^^ J 

du dv {A:'''^ {up, +xj_)Ay^{vp,+y^))^^^ ^ -4as —{x±\^^U'^'^^\y±) (44) 
Jo J -00 Jq Q; P± ~ Pj_ ~ 

Hereafter we use Schwinger's notations (xj_\F{pj_)\yj_) = Jdp e^^P'^~y^^ F{pj_) (the scalar product of the four- 
dimensional vectors in our notations is a; • y = — |(a;*y» + x.y,) + {x,y)±). Note that the interaction with the 
shock wave does not change the a-component of the gluon momentum, same as for the scalar propagator (j3ip 

Formally, the integral over a diverges at the lower limit, but since we integrate over the rapidities Y > Y2 in the 
leading log approximation, we get {AY = Yi — F2) 

g^TduT dv {A:^''^iup,+x^)Ai'''^{vpi+yi^))p,^^ - _ 4a,Ar(.T^l4[/"^4|y^) (45) 

Jo J -00 Pj_ Pj_ 



and therefore 



13 23 



(hereafter = —if "■""'■). The contribution of the diagram in Fig. [Sja is obtained from Eq. pS)) by the replacement 
T"-Uzi ® ^'"^^22 ~^ Uz^T^ <E) Ul^T"', Z2 ^ z\ and the two remaining diagrams are obtained from Eq. |45] by taking 
Z2 — z\ (Fig. int) and z\ ~ z^ (Fig. [SJl). Finally, one has 

(Tr{?7^ C/t^^n)Fig[5] = ^ / rf'^3^[Tr{T"C/,,C/t,r''C/,3C/tj - ^Tr{t/,,C/tj] (47) 

There are also contributions coming from diagrams similar to Fig. [5] but without the gluon-shockwave intersection. 
These diagrams are proportional to the original dipole TrjC/zj^C/j^} ^'^'^ therefore the corresponding term can be 
derived from the contribution of Fig. [S] graphs using the requirement that the r.h.s. of the evolution equation should 
vanish in the absence of the shock wave (when J7 = 1). It is easy to see that this requirement leads to 

(Tr{L^i:^C/irn> = ^^\d^z^^\T,{T'^UzMXU.M^-No^^{U..Ul}\ 

~ J ^13^23 

which gives the BK equation for the evolution of the color dipole in the adjoint representation:. 



d 

dY 



-TriUlui:-} = ^ / d',,^[Tr{T^uluiyT^uluir} N.Tr{Ulul^}] (48) 

~ J -513^23 
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3. Leading order: BFKL evolution of color dipoles 

Next step is the evolution of color dipole. To find the amplitude ((20)) in the leading order (and NLO as well) it 
is sufficient to take into account only the linear evolution of Wilson-line operators which corresponds to taking into 
account only two gluons in the t-channel. The non-linear effects in the evolution (and the production) of t-channel 
gluons enter the four-current amplitude (|20p in the form of so-called "pomeron loops" which start from the NNLO 
BFKL order. With this two-gluon accuracy 

where W^lzi, Z2) is a color dipole in the adjoint representation, see Eq. (|40l) . The BFKL equation for ZY°'(zi, Z2) takes 
the form (recall that a = e^) 

d f 

(J—W{zi,Z2) = (fz3(fz4K-Lo{zi,Z2;z3,Z4)U'^{z3,Z4) (49) 

where 



ifLo(^i, ^2:^3,^4) = ^rfl4!(fH)+fl4!(fHl)_ 5\z,,)5\z,,)ld'z-, ^ (50) 



(Zl - 2)2(2:2 - z)2 



The solution of this equation is easily formulated in terms of lA"{v^za) - projection of color dipole on Lipatov's 
eigenfunctions (|35|) : 

W{v,zo) = (a/ao)"('''z:/^°(i^,Zo) (51) 



4-. LO: amplitude 

The last step is a matrix element of the color dipole operator U°'°{i', zq) "between scalar states", i.e. the correlator 
of color dipole and "bottom pair" of scalar operators 0{x') and 0{y') in the leading order in perturbation theory. 
The easiest way to get this matrix element is to write down the high-energy OPE for the bottom pair of operators 
similar to Eq. ([^ 

Here = pi + + z'^j_, 60 = k'^^ + ie = sOr'^-l'y + ^^'^ 

/i ./ -iio 1020 



where the operator is made from the dipoles V'^{z[,Z2) = 1 — ^ Tr { V^,^ } (cf. Eq. (00])) ordered along the 
straight line || p2 with the rapidity restriction 



V,^=Pexp~ig dup'^A^^iup2 + xj,) , 4(x)= / d^/c 0(A - |/3fe|)e*"^^(fc) (54) 

J — OQ J 

Similarly to the case of the upper impact factor discussed above, the cutoff A for integration in Eq. ([M)) should be 

chosen of order of characteristic /3's in the lower impact factor so ^ ^x''}y')2 ■ 
In the leading order in perturbation theory 

/7y/ \\ "sTT^iV^ (^zi-Wi)'^{z2-W2)'^ 

{Uiz,,Z2)V{wi,W2)) = -—^-2-^ In ) (2 55 

2(7V^^-1) [zi - W2r(Z2 - wiY 

which will be true in the LLA as long as the a and /3 cutoffs do not allow large logarithms In where fc^ is 
characteristic transverse momentum in gluon ladder describing the BFKL evolution. (This is similar to taking /i^ 
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around 1 GeV for the initial point of the DGLAP evolution so the logarithms In ^ can be neglected). Thus, if we 

choose the final point of evolution (|5T|) to be (To ^ ^ {\x — y\±\x' — y'\±XBs)^^ , the correlator of color dipoles 
(W^" {zi,Z2)V^^ (wi, W2)) wiU be given by Eq. ([55]) which translates to 



167r2 



S{zo - z'Q)S{iy + v') + 



I \ 2-4iu 



- 1 J/2(l + 4l/2)2 L - ■ TT\Zo - 

where we used the following orthogonality condition for the eigenf unctions ([35]) . see Ref. [ijj 

^4 



d'^zid'^Z2 j^* , I / \ 

4 E^,^^{zi - Zq, Z2 - ZQ)E^^n{zi - Zo,Z2 - Zq) = 



(56) 



(57) 



6{iy ~ i'')Sm,nS'''^HzQ - Zq) 

+ b{v-rV)bra,-n[.Z^-Z^) [Zo - Z^) I — + j p( 1 + |»| , .^.^p^H .^.^ 

With our choice of o-q for the endpoint of the evolution of the color dipole ([51]) the correlator of two color dipoles 
U'^^ {ly, zo) and (;/', z^) takes the form 



{U''-{,.,zo)V^-{iy',z',)) 



(58) 



iV2 - 1 1/2(1 + 4z/2)2 Vs(a; _ j,)i(x' - y')i 



, ,,,, 2^-^"'5{v-v')T{\+w)T{l-w) 



TT\Zq - ZqI 



r(zz.)r(i 



Combining now Eqs. (|38|) . (|52p . and the above equation we see that the leading-order amplitude is given by Eq. ([23 
with wo(j^) = ^^x{^) and 



i^o(i^) 



iV^ 47r4Q,2 



- 1 cosh^ TTi/ 

which agrees with the leading-order impact factor calculated in Refs. [l^lsij. Here we used the integral 



[(zO-4)2]1-2-L(-2ac'.C^)2 



{-2k' ■ Cor 



(59) 



(60) 



C. NLO: Rapidity factorization 

Now we repeat the same four steps of operator expansion at the NLO accuracy. A general form of the expansion 
of T-product of the currents 0{x) and 0{y) in color dipoles looks as follows: 

T{6{x)d{y)} - Jd^z,d'z2l'^''{z,,Z2)Tv{Ulul^} 

+ 1 d'z,d'z2d'zs /^^°(zi, Z2, zs)[TT{T-Ulul^r^Ulul^} - N,Tr{Ulul^}] (61) 

(structure of the NLO contribution is clear from the topology of diagrams in the shock-wave background, see Fig. [Sj 
below). 

The NLO impact factor for two currents is given by the diagrams shown in Fig. [6] The gluon propagator 
in the shock-wave background at a;* > > in the light-like gauge ^2^;^ = is given by Eq. ([i5)) To calculate 
the next-to-leading impact factor we also need the three-point scalar-scalar-gluon vertex Green function (vertex with 
tails) which is proportional to 



d^z 



— d 

(x — z)2 ^ (z — j/)2Jz^ x'^y'^{x — y) 



^ o J/ = 4i7r " ^ (62) 



Using this formula,, one obtains the contribution of Fig. [nti,b diagrams in the form (the details of the calculation are 
presented in Ref. [IJ]) 



{T{d{x)d\y)}r^^-^ = "U22 I I d'z,4^Tr{T-U.,UlT-U..Ul} 

[N^ -l)TT'^xiyiJ Z(Zi J zf^z^^ Jo a 



e''"323 (63) 
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Let us discuss now the contribution of Fig. [5J:,d diagrams. Since this contribution is proportional to Tr{UziUl^} it 

can be restored from the comparison of Eq. (1551) with the pure perturbative series for the correlator {T{0{x)0{y)}). 
If we switch off the shock wave the contribution of the Fig. |6] diagrams is given by the second term in Eq. ((63)) 
(with U, repaced by 1). On the other hand, perturbative series for the correlator {T{0{x)0{y)}) vanishes [s^ and 
therefore the contribution of the Fig. |BJ:,d diagrams should be equal to the second term in the r.h.s. Eq. ([55)) with 
opposite sign. Thus, we get 



{T{d{x)dHy)}fj^^ 



(iV2 _ l)n^xlyl 



^1 ^2 



d Z3- 



'■13-'23 



Tr{T"C/,,[/t T"C/,3C/t } _ iV,Tr{[/,, (7^ } 



'da 



a 



(64) 



The integral over a in the r.h.s. of Eq. ()64l) diverges. This divergence reflects the fact that the r.h.s. of Eq. (|M)) 
is not exactly the NLO impact factor since we must subtract the matrix element of the leading-order contribution. 
Indeed, the NLO impact factor is a coefficient function defined according to Eq. (|CT)) . To find the NLO impact factor, 
we consider the operator equation ()6ip in the shock-wave background (in the leading order {Uz3)a — Uzs): 



(65) 



{T{d{x)d{y)})A ^ Jd'z,d'z2 l'^''{x,y;z,,Z2){TT{Ulul^}}A 

d^z,d^Z2d^Z3 l'''^''{x,y;z,,Z2,Z3;Y)[TT{T^Uz,UlT^Uz,Ul} - 7V,Tr{C/,, L/tJ] 

The NLO matrix element {T{d{x)d{y)}) a is given by Eq. ^ while 

'd^ zid^ Z2 ix,y;zi, za) (Tr{ Ul Ut^})A 

as follows from Eqs. (15^ and (H5)) . The a integration is cut from above by cr = in accordance with the definition 
of operators (|5D)) . Subtracting (|55)) from Eq. (|M)) we get 



rNLO 



{x,y;zi,Z2,Z3;Y) = 



as{x^y^y 



n\N^ - 1) zf,zi,ZfZl 



'da 



a 



' da 



a 



as{x^y^y 



'13 



In — Zg h C 

4 ^ 2 . 



^4(Ar2 _ 1) zl,zl,ZlZl 
Let us rewrite the operator expansion (1611) in the explicit form [l^ : 

1 



{x-yfT{0{x)d^{y)} = 



^2(7V2 _ 1) 



d zid zo „T 



-12 



(67) 



(68) 



'12 



2 ,2 



Zi riZ- 



13-'23 



ln^aZ3-|+C 



Note that the l.h.s. of the Eq. (|55)) is conformally invariant while the coefficient function in the r.h.s. is not (due to 
[In fcrZa factor). The reason for that is the cutoff in the longitudinal direction (|5(I)) . Indeed, we consider the light-like 
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dipoles (in the pi direction) and impose the cutoff on the maximal a emitted by any gluon from the Wilson lines. 
Formally, the light-like Wilson lines are Mobius invariant. However, the light-like Wilson lines are divergent in the 
longitudinal direction and moreover, it is exactly the evolution equation with respect to this longitudinal cutoff which 
governs the high-energy behavior of amplitudes. At present, it is not known how to find the conformally invariant 
cutoff in the longitudinal direction. When we use the non-invariant cutoff we expect, as usual, the invariance to 
hold in the leading order but be violated in higher orders in perturbation theory. In our calculation we restrict the 
longitudinal momentum of the gluons composing Wilson lines, and with this non-invariant cutoff the NLO evolution 
equation in QCD has extra non-conformal parts not related to the running of coupling constant. Similarly, there 
will be non-conformal parts coming from the longitudinal cutoff of Wilson lines in the J\f = A SYM equation. We 
will demonstrate below that it is possible to construct the "composite conformal dipole operator" (order by order in 
perturbation theory) which mimics the conformal cutoff in the longitudinal direction so the corresponding evolution 
equation has no extra non-conformal parts. This is similar to the construction of the composite renormalized local 
operator in the case when the UV cutoff does not respect the symmetries of the bare operator - in this case the 
symmetry of the UV-regularized operator is preserved order by order in perturbation theory by subtraction of the 
symmetry- restoring counterterms. Following Ref. [31 we choose the conformal composite operator in the form (fTT|) 

J ^13^23 *^13^23 

where a is an arbitrary constant. It is convenient to choose the rapidity-dependent constant a ae~^^ so that the 
[Tr{?7^^C/j^}]™" does not depend on F = Incr and all the rapidity dependence is encoded into a-dependence: 

= Tr{U;,Ut:}+ ^ d'z3^mr'U^,Ul:T-U^^Ul^}-N,Tr{U^^Ul:}]\n ^ + 0{al) (69) 

J ^13^23 ^^13^23 

Using the leading-order evolution equation p8t it is easy to see that -^[^"^{U ziUl^}^"^^ = (with our O(a^) 
accuracy). 

Rewritten in terms of conformal dipoles (|69p . the operator expansion (|68|) takes the form: 



- S/^'^3;^(ln-5^Z|--z^ + 2c)[Tr{r"{7,,t/t,rn;y^^{>t^ (70) 

ZTT J 213^23 ^ ^Zl3'^2?, ' ' 



We need to choose the new "rapidity cutoff" a in such a way that all the energy dependence is included in the matrix 
element (s) of Wilson- line operators so the impact factor should not depend on energy ( = should not scale with p as 



p — > oo). A suitable choice of a is given by oq = k + ie = f^^^'yi + is so we obtain 



/ j2 



o_-2 / " "3 ^2 2 
ZTT ./ ^13-223 



7r2(Ar, 

2 2 

' , -X*y*Zj^2 p-2 , r,^ 



^12 



[TT{rVlUl:T-U:.^^Ul:} - A^,Tr{^,-^;7t;}]} (71) 



where the conformal dipole [Tr{f7j;f7]^^}]™"f is given by Eq. 1^ with uq = 

Now it is evident that the impact factor in the r.h.s. of this equation is Mobius invariant and does not scale with 
p so Eq. (j69p gives conformally invariant operator up to a'j. order. In higher orders, one should expect the correction 
terms with more Wilson lines. 

With the two-gluon accuracy integration over one z in the r.h.s. of Eq. (j71l) can be performed explicitly so that 
the resulting operator expansion takes the form 

{x^y)^T{d{x)dHy)} (72) 
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We need the projection of the T-product in the Lh.s. of this equation onto the conformal eigenfunctions of the 
BFKL equation with spin (of. Eq. 



1 f dzi dzi r, r a^Nr 

Zfn 27r 



In' 7^ - - 2C( ln7^ - ^ + 2) + 2Li2(l - 7^) - y ] } [ 



(-2«-Co)^ 



(73) 



where 



<i>iH = -2^'(i+^.)-2V/(i~^.) + ^ + ^^|^ + 2Cx(.) 



(74) 



and Co = ^ + 2oJ-P2 + zoj^. 

Now, using the decomposition (j37|) of the product of the transverse (5-functions in conformal 3-point functions 
Ei,^n{zio,Z2o) we obtain 



{x-ymd{x)dHy)} 



= — I dv I Szq 



l-2w) V(-2K-Cn)2/ ^ 



47rcosh7ri^ r(l - 2w) V(-2k-Ci 



27r 



<i>iM]4™f('^,^o) 



(75) 



where 



"confl^i^Oj = ^ 



'12 



/ ^2 \ 

(72^) ^conf(^l:^2) 
^ -'10-^20 ^ 



(76) 



is a conformal dipole in the zq, representation. 



D. NLO BK in A/" = 4 SYM 



The typical diagrams for the NLO evolution of color dipole are shown Fig. [T] Here solid line depicts either scalar 
particle or gluino. The contribution of these diagrams is calculated using the gluon propagator in a shock- 
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(VIII) 



FIG. 7: Different types of diagrams for the NLO evolution of color dipole. 

wave background (|43p and usual Yang-Mills vertices. The typical contribution has a form of an integral over Feynman 
parameter which is a fraction of momentum a carried by one of the gluon after splitting described by the three- giuon 
vertex. For example, 



{Tt{U.UI})f,,^ = ^InAF^ du uuj d'zd'z' 



X Tr 



where 



be 



(z — z'y-[uz1^ + uz 



{zi^z2)\]u^[t-' r''''^' 



Z'ij 



[z — z')2[mz|3 -|- uz. 



)]}UI (77) 



2 2 2 J , ^ J 

(^13 ~ ^14)9*"' + —Zd.i^U + ~^13^34' 



2 

u 



71^3 



(2 2 \ ij : ^ i 3 I i 3 

1^23 ^24^3 + ^34^24+ -^23^34 



(78) 



and u = 1 — u. It is easy to see that integral ([77]) diverges as u and m — > 1. This divergence comes in the 
momentum space from the integrals of the type 

daida2 — tt^ -TTl T 

where aipi and ki± are longitudinal and transverse momenta carried by gluons in the loop in Fig. [7j. If we put a 
lower cutoff > a' on the integrals we would get a contribution ~ In^ -p- coming from the region 02 3> ai > cr' (or 
ai ^ a2 > (j' ) which corresponds to the the square of the leading-order BK kernel rather than to the NLO kernel. To 
get the NLO kernel we need to subtract this (LO)^ contribution. Indeed, the operator form of the evolution equation 
for the color dipole up to the next-to-leading order looks like 



d 



dY 



(79) 



where Y = Incr. Our goal is to find -J^nlo by considering the l.h.s. of this equation in the external shock-wave 
background so 



Shockwave — Shockwave 



ckwave 



The subtraction (f5(I)) leads to the ^ prescription (p^ for the terms divergent as ^ (similarly, 4 



(80) 
for the 



contribution divergent as m — ^ 1). With this prescription, the integrals over the Feynman parameter converge. The 
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calculation of diagrams is carried out in Refs. [13, [l^ and the result has the form 



dY 



d z-i 



2 r 

[Zoo ^ 



'13^23 



47r 



2 In ^ In 



13 i„ -'23 



-\2 ^12 



4r^ 



2 '' d'^Z^d'^Zi 212^34 



"34 



"13^24 

X 



'12^34 



-13^24 ^23-^14 



In 



'13^24 



'TY\bb' 
Z3I 



(81) 



All terms in the r.h.s. of this equation are Mobius invariant except the double-log term proportional to In ^ In 
As was discussed in the Introduction, the reason for this non-invariance is the cutoff in the longitudinal direction 
which violates the formal invariance of the non-cut Wilson lines. 

It is worth noting that conformal and non-conformal terms come from graphs with different topology: the conformal 
terms come from 1— >3 dipoles diagrams (see Fig. 6 in Ref. [l^) which describe the dipole creation while the non- 
conformal double-log term comes from the 1—^2 dipole transitions (Fig. 9 in Ref. [l^)) which can be regarded as a 
combination of dipole creation and dipole recombination. 

Our aim is to rewrite the evolution equation in terms of composite conformal dipoles (|69p . In the next-to-leading 
order the conformal dipole has the form 



(82) 



as 
27r2 

„2 



d z-i 



'-12 



'13^23 



47r 3 



[Tr{T^UluZT^UluZ} - N,Tr{UluZ}] In 



4az 



12 



^2 ^2 



"s / t2_ t2^ ^12 foi„ ^12^34 I [1 I ^12^34 1 i fisf^ 1 i 

d Z3d Z4-^-^—^<2ln-^—^ + 1 + ————— in-2--2-|ln--2-/(z,) 
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where function / is undetermined in the NLO (to fix it one needs the NNLO accuracy). Using the equations 
for ^Ti{UY,UZ} and ^[Tr{T''UlUl^T''U,.^Ul^} - 7VeTr{/7]^[/]f }] from Ref. [H one can demonstrate that 

^[Tr{;7.,;7tj]:"^ =^0. 

Differentiating now with respect to a we get 
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where 
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^14^24 '''-^ ^14''24 



is a conformal y-independent operator found in Ref. [14] . One sees now that the evolution equation with respect to 
parameter a (j83p is obviously Mobius invariant. 



E. NLO evolution of the conformal dipole 

With out two-gluon accuracy the evolution equation reduces to 



2a— ZY°o„f(zi,Z2) = J d zsd Z4[Ki^o{zi,Z2;z3,Z4) + Kjqi^o{zi,Z2;z3,Z4)]U^^^f{z3,Z4) (85) 
where the kernel K-^^i^oizi, Z2', 23, Z4) in the first two orders has the form [l^, [23] (see also Ref. [3 
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Once we know that the kernel -RTnlo is conformal, its eigenfunctions are fixed by conformal invariance (see Eq. psp) 
and calculation of the eigenvalues reproduces the pomeron intercept (l?f)) . 

d^Z3d^Z4 K{zi,Z2\Z3,Z4)Ey^n{Z'in,Z4n) = u)[n, v)Ei,^n{zw, Z20) (87) 

For the composite operators with definite conformal spin (j76p the evolution equation (|85p takes the simple form 

2a^Zi,Vt(^,--o) - c.(j.)Z:/,Vf(^,--o) (88) 

The result of the evolution (|88p is 

^c"o„f(i^,^o) = {araY^-^''^KUi^,z^) (89) 

where the endpoint of the evolution a should be taken from the requirement that the amplitude of scattering of 
conformal dipoles with "normalization points" a and b should not contain large logarithms of energy so it will serve 
as the initial point of the evolution (cf. Eq. ([55)) for the leading order). 

F. NLO scattering of conformal dipoles and the NLO amplitude 

The last step of our program is the "matrix element" of the conformal dipole 

{T{UtA^,z,,)0{x')0\y')}) (90) 
As we have done for the leading order, we start with the high-energy OPE for 0{x')0'^ [y') 
{x' -yrT{d{x')&{y')) 
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where 60 = l/re'^and the conformal operator V^onf (^i' -^2) is given by Eq. (l53l) . 

Now we multiply Eq. ([?T|) by the NLO amplitude of scattering of two conformal dipoles [s^ : 
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From this equation we see that we need to stop the evolution ((5^ at a = 1/h. With this choice of a the correlation 
function takes the form: 
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(93) 
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Finally, substituting this amplitude in Eq. (1751) we obtain Eq, 
Nl 



1 



27r 



iV2 - 1 cosh^ 



1 



2tt 



with 
x(^)(4C + 



1 + 4i/2 



4C 



47r4a2 



- 1 cosh^ -. 



1 



a,Nr_ 



o(«?)} 

0(aD} 



(94) 



which gives the pomeron residue in the next-to-leading order. 

It is instructive to represent this result in a way symmetric between projectile and target as a product of spectator 
impact factor, target impact factor, and scattering of two (conformal) dipoles as shown in Fig. [51 
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FIG. 8: High-energy factorization into the product of two impact factors and dipole-dipole scattering 

Combining the operator expansions (|75p and (I?!]) we get 
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The scattering of two conformal dipoles is obtained from Eq. 
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Substituting this amplitude with gq = k ^ + ie and bo = k' +ie in Eq. (1951) we reobtain Eq. (1231) with the pomeron 
residue (|M1) . 

V. NEXT-TO-LEADING ORDER IN QCD 

In A/" = 4 SYM we were able to perform all four steps of our program and arrive at explicit result (l^5|) with uj{i') 
and F^f) given by Eqs. ([3]) and (j94| . respectively. In QCD, we do not have conformal invariance so there is no 
general formula for the pomeron contribution similar to Eq. (|23p . Moreover, by the same reason (absence of Mobius 
invariance) we cannot solve the equation for evolution of color dipoles analytically. Here I will describe the first two 
steps of OPE program - (1) operator expansion in color dipoles and (2) evolution equation for color dipoles. 



A. High-energy OPE and photon impact factor 



I will consider the process of deep inelastic scattering (DIS) at small Bjorken x. As I discussed in Sect. Ill, the 
virtual photon splits in to qq pair which moves in an "external" field of target's gluons. To calculate the leading-order 
impact factor, we need the quark propagator in the shock-wave background. It has the form 0] 



27r2(a;-z)4 '27r2(j/-z)4 



(97) 





-a^ 

factor is a product of two propagators ((57|). see Fig. |3]with solid lines representing quarks 



where ij: = —iip'^jQ. We use Dirac matrices 7^^ = — ^ ) where cr'^ are defined in Eq. (fT8|). The LO impact 
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Hereafter I use tr{...} to denote the color trace in the fundamental representation (and reserve the notation Tr{...} 
for the trace in the adjoint representation). The above equation is explicitly Mobius invariant. In addition, it is easy 
to check that T5^(r.h.s)=0. 

ox ^ ^ ' 

Feynman diagrams for the photon impact factor are the same as in Fig. [5] but with the solid lines representing 
quarks. The calculation of the NLO impact factor is similar to the scalar case. The master formula is (cf. Eq. (|62|) ) 
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which gives the 3-point ipipF^^, Green function in the leading order in g. Using this formula one easily obtains 
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Performing integrals over z, 's and taking traces one gets after some algebra the NLO contribution in the form (recall 
that z2. = -2(C, • Cj) and Z, = -^{k- 
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where we have promoted the equation to the operator form. It can be demonstrated that g|jr[r.h.s. of Eq. (jlOip ] = 
which reflects the electromagnetic gauge invariance. 

As for the N=4 case, it is convenient to re-expand r{jp(a;)j,y(y)} in composite operators obtained by replacement 
^ = ^ (and Tr tr) in Eq. ^ 
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where dots stand for terms not displayed here for brevity. In QCD, this composite operator is no longer conformal 
but still more convenient than the original dipole since the evolution equation for this operator splits into the sum of 
the conformal part and the running-coupling part proportional to 6 = ■^Nc — |n / (see the next Section). The final 
expansion of T{j^{x)j^{y)} in composite dipoles (|102p has the form (cf. Eq. ([71]) ) 
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(103) 



where we set uq — k 



ie same as for the A/" = 4 case. The explicit expression for four- Wilson-line 



composite operator [tT:{Uz-LUl^}ir{tJz3tjl^}\a can be obtained from Eq. (j84l) by usual substitution T" — >• i° and Tr — >• 

tr. (It is worth noting that at large \^'':{UziUl^}ir{lJz^tjl^}]a — [^T:{UziUt^}]a[^T:{Uz3Ul^}]a- ) It is easy to see now 
that the NLO impact factor ()103p is Mobius invariant. 
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B. Evolution of color dipoles in QCD 



Second step is the evolution equation for composite operator (|102[) . The types of Feynman diagrams (in the shock- 
wave background) in QCD are the same as in A/" = 4 SYM (see Fig. [7]) but now the soHd hues in Figs. [7] VII and 
VIII denote quarks rather than scalar or gluinos. The result of the evolution of the color dipole with rapidity cutoff 
(El is [H: 
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where we use the MS scheme. The NLO kernel is a sum of the running-coupling part (proportional to b), the non- 
conformal double-log term ~ In In (same as in A/" = 4 case) and the three conformal terms which depend on 

■^13 "^13 

2 2 2 2 

the two four-point conformal ratios '^l^^?/ and . 

A natural guess for QCD result as opposed to A/" = 4 answer would be some sort of tree-level conformal structure 
which "get dressed" by the running coupling constant. At the NLO level, this would correspond to the sum of the 
conformal part and the running-coupling part proportional to b. As we see, the Eq. (jl04p does not quite look like 
such sum due to the presence of the double-log term. It turns out, however, that if one rewrites the evolution equation 
(|104l) in terms of composite operators (|102p . the NLO kernel separates into the sum of conformal and running-coupling 
parts [13] : 
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Following the analysis of Ref. [l^ I will outline how the above kernel reproduces the NLO BFKL eigenvalues Q (for 
details see the Appendix). 

In the two-gluon approximation we get 



(106) 



24 



where 
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is the color dipole in the fundamental representation. In this approximation the composite dipole (jl02p reduces to 
(cf. Eq. dMl)) 
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It is easy to see that ^ ( = (t-^ ) of the r.h.s. vanishes. 
The evolution equation (jlOSt turns into 



2a—lla{zi,Z2) 

da 



27r2 



d'^3 



'13^23 



47r 
a'iV; 



2 

s ^ ' c 



'ln£i2M^ 
^ 4 
;2. 



2C) 



'23 



'12 



1 ^13 ^67 



10 



\Ua{zi,Z3) +Ua{z2,Z3) -Ua{zi,Z2)] 



fl + 



d Z3d Z4 [^ 2:12^34 ^12^34 



"34 



'13^24 



-14^23 



'12^34 



-13^24 



12^34 



14^23 



In 



^13^24 



3^12^34 



Zi ^2:9^ Zi ^ Zo 



, ^13^24 
2 2 
^14^23 



-y2 -y2 I -,2 -,2 _^2 ^2 ^2 _^2 

^13-'24 "T ^14-'23 ■'12^34 ^13^24 

.-^2 ^2 '^2 

^13^24 -'l4-'23 -'l4-'23 



where we used formula [14] 
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For the case of forward scattering {U{x,y)) — U{x ~ y) and the linearized equation p09p can be reduced to an 
integral equation with respect to one variable z = zyi- Using integrals (104)-(106) from Ref. 19] and the integral 
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FIG. 9: Renormalon bubble chain of quark loops. 
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The function — (^z^rp-ln^ ^ + F{q,q') + $(q, g') enters the NLO BFKL equation in the momentum space 0] and 
since the eigenfunctions of the forward BFKL equation are powers both in the coordinate and momentum space, it is 
clear that the corresponding eigenvalues coincide. Moreover, it can be demonstrated explicitly that 



1 



— \n'^+F{q,q') + ^q,q')] = - ^ In^ ^ + F{z, z') + ^z, z') (114) 
} q [z z ) z 



{q-q' 

so the conformal (A/" = 4) part of the forward kernel looks the same in the coordinate and in the momentum 
representation. As to the running-coupling part (the first term in the r.h.s. of Eq. (|120p ). in the Appendix we 
demonstrate that it also agrees with the eigenvalues ([3]). 



C. Argument of the coupling constant in the BK equation 

I will not discuss here the steps 3 and 4 of our OPE program. The reason is that in DIS from nucleon or nucleus 
the evolution of color dipoles is non-linear and the analytic solution is not known at the present time. Even in the 
simpler case of forward 7*7* or onium-onium scattering where the dipole evolution is described by linear NLO BFKL 
equation, it is impossible to solve this equation since we do not know the argument of the coupling constant. Moreover, 
it is not known how to solve analytically this equation even if we take some simple model for the argument of coupling 
constant like the size of the parent dipole. 

Still, the first step towards the solution would be to figure out the argument of coupling constant in the NLO 
BFKL equation. To get an argument of coupling constant we can use the renormalon-based approach (for a review, 

see Ref. :3 73) and trace the quark part of the /3- function proportional to n/. In the leading log approximation 

2 

In ^ '-^ 1, tts ^ 1 the quark part of the /3- function comes from the bubble chain of quark loops in the shock- wave 
background. We can either have no intersection of quark loop with the shock wave (see Fig. [9^) or we may have one 
of the loops in the shock- wave background (see Fig. Wp)- 

The sum of these diagrams yields 

X ^(l-^ln^)diU^\l) 1 ^ (115) 

p2(i + |.ln^)V 6n ^ + 
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where we have left only the /3- function part of the quark loop. Replacing the quark part of the ^9- function —^Uf In ^ 

2 

by the total contribution j^bln^ we get 

(p — t)^ 

In principle, one should also include the "renormalon dressing" of in Eq. (|105p . We think, however, that they 
form a separate contribution which has nothing to do with the argument of the BK equation. 

To go to the coordinate space, we expand the coupling constants in Eq. (jll6p in powers of ag = agip^), i.e. return 
back to Eq. (jllSp with f^rif — )■ —bj^. Unfortunately, the Fourier transformation to the coordinate space can be 
performed explicitly only for a couple of first terms of the expansion as(p^) ~ as — In p^/fi^ + lnp^//i^)^. 
With this accuracy [H, HI] 

±Tr{U^^Ul} = I d'z [Tt{U,,uI}Tt{U^,uI} - N,Tr{U^,Ul}] (117) 



(118) 



-^13^23 ^13 ^0!s(223) ' ^23 ^'-'^s(^13) 

where dots stand for the remaining terms. (Here we promoted Wilson lines in the r.h.s. to operators). 
When the sizes of the dipoles are very different the kernel of the above equation reduces to 



1^:121 < W-i\ 



27r 2:^32:23 

kl3«|zi2|,k23| 



^^i^ k23|«Nl2|,Nl3| (119) 

so the argument of the coupling constant is the size of smallest dipole. The numerical approach to solution of the the 
NLO BK equation with this running coupling constant is presented in Ref. (39j . 



VI. CONCLUSIONS 



The main conclusion is that the rapidity factorization and high-energy operator expansion in color dipoles works at 
the NLO level. There are many examples of the factorization which are fine at the leading order but fail at the NLO 
level. I believe that the high-energy OPE has the same status as usual light-cone expansion in light-ray operators so 
one can calculate the high-energy amplitudes level by level in perturbation theory. As an outlook my collaborator G. 
Chirilli and I intend to apply the NLO high-energy operator expansion for the description of QCD amplitudes. (The 
intermediate result for the impact factor (|103|) is the first step of that program) . There are many papers devoted 
to analysis of the high-energy amplitudes in QCD at the NLO level (see e.g. Refs. [1, H^, |4l|) but all of them 
use traditional calculation of Feynman diagrams in momentum space. In our opinion, the high-energy OPE in color 
dipoles is technically more simple and gives us an opportunity to use an approximate tree-level conformal invariance in 
QCD. Moreover, the exact prescription for separating the coefficient functions (impact factors) and matrix elements is 
somewhat tricky in the traditional approach while it comes naturally in the framework of OPE logic. When finished, 
the calculation of the photon impact factor for the structure function ^2(2;) of deep inelastic scattering will compete 
the analysis of the small- a; behavior of DIS structure functions at the NLO level. The study is in progress. 

The author is grateful to L.N. Lipatov for more than thirty years of valuable discussions and guidance. This work 
was supported by contract DE-AC05-06OR23177 under which the Jefferson Science Associates, LLC operate the 
Thomas Jefferson National Accelerator Facility. 
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VII. APPENDIX: COMPARISON TO NLO BFKL IN QCD 



We will compare our evolution equation for to similar equation obtained from the NLO BFKL momentum-space 
analysis [3, The linearized version of the Eq. (I104p has the form (cf Eq. 



asNc 
27r2 



{z-z'f-z''^ , {z^z'Y . .67 7r2 



In- 



(z — z'Y'z 

+ ( 



z^a^ 



10 



z 



aiN] 



47^3 



9 3'''^ 9 
d'z' ^ [F{z, z') + <i>(z, z')] (z') + ^C(3)Z^^ (z) 



[W{z~z')+W{z')-W{z)] 



where F{z,z') and $(z,z') are given by Eqs. (|112l) and (I113p . respectively. 

To compare the eigenvalues of the Eq. (I120p with NLO BFKL we expand {x, 0) in eigenfunctions ([35 



(ZY^(x^,0))= ^ 



— -ir +'iOO 



n— — oo 2 



— —zoo 



27ri 



(120) 



(121) 



compute the evolution of {U{n, 7)) from Eq. (|120p and compare it to the calculation based on the NLO BFKL results 
from [1, [25 . (Here we will not consider the quark part of the NLO BK kernel - the agreement of that part with the 
Hf term in the NLO BFKL kernel was proved in Ref. [42*1). 
The relevant integrals have the form 



i-|d2z[2(zVx^)V"*-l](^^^ - Xinn) 



<fz [2(zVx2)Te"* - 1] 



— z)2z' 
1 1 

(a; — z)2 z2 



In 



(a; — z)2 



d^z (z7a;2) 



In 



(x — z)2 , z2 



(.T — z)2z2 ^ 

where x(n, 7) -2C - -0(7 + f ) - -0(1 - 7 + f ), and 
i /d^z' (z'V^')^^'e™^'F(z,z') 



■In- 



= oX"(",7) +X'(",7)x(»^,7) 



47x(7) 



(122) 



71/ \ 2 + 377 
iVfJ(3-27)(l + 27) 



77 



"}(1-^ 



cos 7r7 



iV3 / 2(3 - 27)(1 + 27) J (1 - 27) sin2 7r7 



F(n,7) 



where [28 



I d^z {z'"^ /z'^y-^e""^' ^{z,z') = -$(n,7)-$(n,l-7) 
27r ./ 



00 / 00 



(123) 



(124) 



(note that x(0, ^ + iv) = x{^) and $(0, ^ + iu) = <i>(zy), see Eq. Q). The convenient way to calculate the integrals 
over angle (j) is to represent cos n(j) as r„(cos(/)) and use formulas for the integration of Chebyshev polynomials from 
Ref. [28)]. 

Using integrals (|122p - (I123P one easily obtains the evolution equation for U{n,j) in the form 



dY 



(Zi^(n,7)) 



TT I 



l-^(A+in4-2C) + ^^^f^-^-HZ^) 
47r M7 47r V 9 3 9 J 



27X('^,7) 



asN_c r 
47r 



X" (n, 7) - 2x(n, 7)x'(«, 7) + 6C(3) + F{n, 7) - 2$(n, 7) - 2$(n, 1 - 7)] } (^"^ 7)) 



(125) 
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where x'(n,7) = ^x("-,7) etc. 

Next we calculate the same thing using NLO BFKL results [2^ . The impact factor (<?) for the color dipole 
U[x,y) is proportional to as{q){e'^'^^ — e^'^y){e^^'^^ — e~"^^) so one obtains the cross section of the scattering of color 
dipole in the form 



1 f d^q d'^q 



47r^ J g2 



a-\-ioo 



as{q){e^'^^ - l)(e-'^- - l)<i>s(g') f ^(—X G^{q,q') 

2,111 \qq' J 



2TTi V qq' 



(126) 



where Gi^{q,q') is the partial wave of the forward reggeized gluon scattering amplitude satisfying the equation 

OjGUq,q')^S^'\q-q')+ [ d'pK{q,p)GUp,q') (127) 



and $_b(<z') is the target impact factor. The kernel K{q,p) is symmetric with respect to q o p and the eigenvalues 
are 



7-1 



as{q) 



X{n,l) + ^^^'5(n,7) 



(128) 



5{n,l) = -^[x'(n,7) + x'(",7)] + (y-y- Y^)x(",7)+6C(3) 
-X"(n, 7) + F(n, 7) ~ 2$(n, 7) - 2$(n, 1 - 7)} 



The corresponding expression for {U{n,^)) takes the form 
1 nn r(-7+f) fd^qd'^q' 



{U{n,^)) 



27r2 ™" 2 r(l + 7+f)y (72 ^,2 
where 6* is the angle between q and x axis. Using Eq. (I127P we obtain 

ds 



-(^) G.(g,g') (129) 



27ri V gq' 



(130) 



1 7m r(-7+^) fd^qd^q' „„. ../g^X-^ , „ r+*°° do; / s \" /" , , . . 



The integration over q can be performed using 

d\aM{^y''e"^'^K{q,p) 



' , has , . , asNg ; 



(131) 



(recall that K{q,p) = K{p, q) and as{p) — as — |^ In ^ with our accuracy). The result is 



r(-7+f) fd^pd^' _ 



s^mnn)) - -^cos.^ r(l + 7 + f)7 p2 ^,2 



\pq' / TT I 2 An 2 An 2 . 



(132) 



;oo 27ri \pq 

where the angle ip corresponds to p. Since w ^ we can neglect terms ^ w in the argument of S and expand 

xin, 7 — ^) ~ x(n, 7) — ^x'(": 7)- Using again Eq. (I127p in the leading order we can replace extra oj by ^Ncx{n, 7) 
and obtain 



.-(^(n,7)) 



as Trn r(-7+f) f d^pd^q' _ / p^ 

fOft / —P ^1 



27r 



X $b((7') 



27ri Vpg' 



2 r(l + 7+f)7 p2 ^,2 vv 

(A)"G.(p,g')^A^c[x(",7) - ^X'K7) + ^[-^(",7) - 2x(n,7)x'(",7)]l (133) 
V vq' / TT I An An J 



Finally, expanding al{p) ~ as{p){as - ^ In ^)as{^l) we obtain 
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s-{U{nn)) = 
as 



asN, Trn r(-7+§) 
■ cos 



27r3 



2 r(i + 7-, 2. 



i^YG^ip^q') (134) 



which can be rewritten as an evolution equation 



An 
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^2 _ "1 
' 4 



2C-ln4- 



a N ~| 



ba, 

asNc 
An 



X"K7) - 2xK7)x'('^,7) +6C(3) + i^(n,7) - 2$(n,7) - 2$(n, 1 - 7)] }(Z^(n,7)) (135) 



47X(7) 



We see that this eigenvalue coincides with Eq. (|125p . 
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